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We study quantum Hall ferromagnets in the presence of a random electrostatic impurity potential, 
within the framework of a classical non-linear sigma model. We discuss the behaviour of the system 
using a heuristic picture for the competition between exchange and screening, and test our conclu- 
sions with extensive numerical simulations. We obtain a phase diagram for the system as a function 
of disorder strength, A, and deviation, Sis, of the average Landau level filling factor from unity. 
Screening of an impurity potential requires distortions of the spin configuration. In the absence of 
Zeeman coupling there is a disorder- driven, zero-temperature phase transition from a ferromagnet 
at small A and \Su\ to a spin glass at larger A or \Su\. We characterise the spin glass phase in terms 
of its magnetic and charge response. 
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I. INTRODUCTION 

Disorder and interactions have competing conse- 
quences in quantum Hall ferromagnets (QHFM). In this 
paper, we study how a disordered impurity potential can 
give rise to a spin-disordered ground state. We also dis- 
cuss the influence of disorder on the magnetic and charge 
response of such a system. We use a classical spin model 
throughout to describe the quantum Hall ferromagnet. 

Coulomb interactions lead to spin correlations in a 
quantum Hall system. For electrons fully occupying the 
lowest Landau level (filling fraction v = 1), exchange is 
responsible for a spin-polarised ground state, even in the 
absence of Zeeman energy. This is the consequence of 
Hund's rule as applied to a macroscopically large num- 
ber of degenerate Landau-level orbitals. The resulting 
quantum Hall ferromagnet is especially interesting as a 
system in which the spin, configuration and the charge 
density are closely linkedEl. At v = 1 and if Zeeman en- 
ergy is large, charge enters the spin-polarised system as 
minority-spin electrons. However, if Zeeman energy is 
small or vanishing, the charged excitation of lowest en- 
ergy is not a bare spin-half electron, but a bound state of 
an electron with many spin waves. In classical terms, this 
occurs because the minority spin polarises its local fer- 
romagnetic background, and the composite object may 
be viewed as a topological excitation, or texture in an 
ordered ferromagnet — a skyrmionH. Similarly, an anti- 
skyrmion, with topological charge of the opposite sign, is 
produced when charge is removed from a filled Landau 
level. In this description, the deviation of local charge 
density from that of a filled and ferromagnetically po- 
larised .Landau level is proportional to the topological 
densityu of the spin configuration. 

In a clean system with sufficiently small Zeeman en- 
ergy, skyrmions or anti-skyrmions can be created at zero 
temperature on varying the average filling factor from 
v = 1 to larger or smaller values. For a disordered quan- 
tum Hall ferromagnet, the coupling of an electrostatic 
impurity potential to the charge density offers an addi- 



tional mechanism by which spin textures arc nucleated 
at zero temperature. The physical consequences of this 
coupling are the subject of this paper, which provides a 
detailed account of work presented in outline in Ref. |[ 

The interplay between disorder and exchange in quan- 
tum Hall ferromagnets has been examined previously 
from several different viewpoints. Fogler and Shklovskiitl 
have developed a theory particularly applicable in-.higher 
Landau levels. Building on earlier discussionso, they 
set out a mean-field treatment in the spirit of Stoner 
theory. For odd integer filling in the absence of Zee- 
man coupling, they find a transition between ferromag- 
netic and paramagnetic ground states with increasing 
disorder strength. They suggest that this transition 
should be apparent in transport measurements, in which 
the ferromagnetic phase is characterised by spin-resolved 
Shubnikov-de Haas oscillations, and the paramagnet by 
spin-unresolved oscillations. Experimentally, a transition 
of this kigd is observed with decreasing magnetic field 
strengtbuQ, and its sharpness suggests that its origin is 
indeed cooperative. 

Within the Fogler- Shklovskii approach, local moments 
are all collinear in the ferromagnet and vanish at the tran- 
sition to the paramagnet. An alternative scenario may 
arise in the lowest Landau level near v = 1, in which the 
QHFM responds to disorder mainly through the direc- 
tion rather than the magnitude of its local magnetisation. 
Some indications that this can happen come from calcula- 
tions for the fully-polarised ferromagnet at weak disorder. 
Here, a reduction in spin stiffness with increasing disorder 
strength has been interpreted by Greenu as a precursor of 
a non-collinear phase. Moreover, even weak disorder may 
nucleate a dilute glass of skyrmions and anti-skyrmions 
at the maxima and minima of the disopdered potential, 
as discussed by Nederveen and NazarovEj and examined 
further in the preseat paper. In addition, Sinova, Mac- 
Donald and GirvinEll have shown that, at intermediate 
disorder strength, both reduced and non-collinear local 
moments emerge from a numerical solution of Hartree- 
Fock theory for a model with Coulomb interactions and 
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spatially uncorrelated disorder, while transport proper- 
ties within a Hartree-Fock theory have been discussed by 
MurthjM 

In general, we believe that the relative importance of 
local moment reduction versus the formation of spin tex- 
tures for dirty QHFM's will depend on the nature of the 
disorder. In this work, we concentrate on textures, which 
are favoured by a smoothly varying impurity potential at 
v = 1. We find that the ferromagnet gives way to a spin 
glass at strong disorder. 

The plan of this paper is as follows. In the next section, 
we describe a spin model with quenched disorder which is 
intended to capture the physics of the disordered QHFM. 
This is followed in section III by a discussion of the phase 
diagram for the model as a function of filling factor and 
disorder strength, using heuristic arguments and scaling 
ideas. The conclusions are supported by the results of a 
Monte Carlo simulation of a lattice version of the spin 
model. The technique is outlined in section [V and the 
resu lts are presented in section |v|. In sections VI and 
VII, we discuss the compressibility and conductivity of 
the system in order to characterise the charge response 
of spin-disordered ground states. 



II. SPIN MODEL 

Consider a two-dimensional electron gas in a strong 
perpendicular magnetic field B, with Landau- level filling 
v close to unity. The electrons are subject to an impu- 
rity potential V(r) and an electron-electron interaction 
energy U(r). As a first step, let us omit the exchange 
interactions and the Zeeman energy. Then the electron 
density p(r) is determined by the balance between dis- 
order and interactions, or in other words, screening. We 
treat this using Thomas-Fermi theacv. Such an approxi- 
mation has been applied by EfrostJIi-a to the comparable 
problem in spin-polarised Landau levels when v lies near 
half-integer values. The ground-state charge density p(r) 
at weak disorder is determined by the condition that the 
Hartree potential should match the chemical potential p 
everywhere: 



p = V(r) + U(r- r>(r')dV . 



(1) 



This approach is valid in the case where the resulting 
local filling factor varies smoothly on the scale of the 
magnetic length, Ib = (h/eB) 1 / 2 , and only has small 
fractional deviations from v = 1. so that 



6p(r) ee p(r) - « p(r) 



(2) 



However, the Thomas-Fermi picture of good local 
screening does not take in account exchange interac- 
tions. Provided electron density fluctuations are small 
and smoothly varying, ferromagnetic exchange leads lo- 
cally to a maximal ferromagnetic polarisation of the 
electron spins. This local magnetisation may vary in 



space. Denoting its direction by the three-component 
unit vector S(r), spatial fmctua^pps in spin orientation 
are linked to electron density by 
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Sp(r) = — eij S ■ [diS x djS 



(3) 



This direct connection is specific to the quantum Hall 
ferromagnet — a varying electron density implies a vari- 
ation in the direction of the local magnetisation and vice 
versa. Such spin textures cost exchange energy. So, a 
proper description of the system must include exchange, 
impurity and Hartree contributions to the total energy 
of a dirty quantum Hall ferromagnet. This brings us to 
study the Hamiltonian 



H = 



J 



\VS(r)\ 2 + [V(r)-p]Sp(v) 



where 



J 



1 



d 2 r, 



(4) 



(5) 



16(2tt) 1 /2 Aneoerh 

is the exchange couplingS. (e r is the relative permittivity 
in the semiconductor.) At this point we have chosen for 
simplicity a short-range Hartree interaction, U(r — r') = 
Uo6(r — r'). We have also absorbed the constant Uq/2ttI b 
into the chemical potential p. 

As mentioned above, we will work with a disordered 
potential V(r) that is smooth on the scale of the mag- 
netic length Ib- For simplicity, our discussion of this 
continuum model will assume a Gaussian distribution 
with standard deviation A and correlation length A much 
larger than the magnetic length I b ■ (Our numerical study 
will use a lattice model with a bounded distribution with 
uncorrelated disorder.) 

In restricting our study to this model, we neglect quan- 
tum fluctuations of S(r). This semiclassical approxima- 
tion is justified for smooth variations, with |V5(r)| <C 
l B \ Our aim in the following is to understand the 
zero-temperature phase diagram of the model defined by 
Eq. (ji|), as a function of disorder strength A and average 
charge density (Sp), the spatial average of 6p(r). We will 
characterise its ground states via their response functions 
and excitations. 

We finish this section by comparing this model with 
some other examples of disordered systems. As an elec- 
tron system, it is unusual in that there is an exchange gap 
for single-particle excitations, even if the ground-state 
spin configuration S(r) does not have long-range ferro- 
magnetic order. This means that the only low-energy ex- 
citations involve collective spin modes. As a ferromagnet 
with quenched disorder, the system is also unusual in sev- 
eral ways. First, the link between spin and charge means 
that the spin system responds to applied electric fields. 
We calculate in the following the wavevector-dependent 
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dielectric constant, e(g), and compare it with behaviour 
found in more conventional disordered electron systems. 
Second, due to the same coupling, spin waves generate 
an electric dipole moment. This means that spin waves 
contribute to the optical conductivity a(u>). More gen- 
erally, the coupling to disorder in this model leaves spin- 
rotational symmetry intact but breaks time-reversal sym- 
metry. This is in contrast to the effect of random Zeeman 
fields, which break both symmetries, and to random ex- 
change interactions, which preserve both symmetries. 

III. PHASE DIAGRAM 

We begin with a qualitative discussion of the phase 
diagram as a function of disorder strength and average 
charge density. We employ scaling arguments to obtain 
the phase boundary for breakdown of long-range ferro- 
magnetic order. 

The Hamiltonian of Eq. (|J) for the continuum model is 
minimised by the spin configuration which satisfies (ap- 
pendix |A|) : 

JW 2 SxS = eijdiVsi^djS (6) 
Vk(r) = V(v)-p+U 8p(r) (7) 

where Vr the local Hartree potential and 8p is defined 
in terms of the local spin configuration by Eq. (||). This 
equation for the ground state configuration shows the 
interplay between exchange and disorder in the model. 
Since it is difficult to tackle the non-linear equation di- 
rectly, we proceed using heuristic arguments instead. 

The model is characterised by two energy scales: the 
exchange energy J and the disorder strength A. There 
are also two length scales: the correlation length A of the 
disordered potential and 

L H = (Uo/ J) 1/2 , (8) 

which we call the Hartree length. The significance of 
the Hartree length can be made clear by considering a 
skyrmion of fixed shape and radius R in a clean system. 
The contribution to its total energy from exchange is ~ J 
and independent of size, while that from Hartree inter- 
actions is size dependent, being ~ Uq/R 2 . Comparing 
these contributions, one sees that exchange dominates 
on length scales large compared with Ljj, while Hartree 
interactions dominate at smaller distances. 

Our central hypothesis is that the competition between 
interaction and disorder in this system is characterised by 
the Hartree length only. In the following, we also use the 
limit Lh > A as a source of simplifications. 

A. Filled Landau level 

Let us consider first the effect of disorder on a system 
at v = 1, imposed by setting (Sp) — 0. Without impu- 
rities (V(r) = 0), the system is a perfect ferromagnet. 



Moreover, there is a thresholdBlBO 

|V(r)| = An J. (9) 

below which an impurity potential is unscreened. It 
arises because, for any S(r), one has |V5(r)| 2 > 87r|<5p(r)| 
and hence 

H> J [AnJ\8p{r)\+V{r)8p{r)]d 2 r (10) 

Thus, if |V(r)| is below the threshold Air J everywhere, 
the ground state is the perfectly aligned ferromagnet with 
Sp(r) = 0. 

At weak disorder (A < J), |V(r)| in most parts of the 
system lies below the threshold. With an unbounded po- 
tential distribution, the ground state spin configuration 
therefore consists of a dilute glass of skyrmions and anti- 
skyrmions, nucleated at rare positions where_fJ^(r)| is 
large, as discussed by Nederveen and NazarovE£l. Away 
from these positions we expect that the ferromagnetic 
order is essentially unaffected by the impurity potential. 
A careful treatment of this regime is however quite sub- 
tle, since the spin deviation due to an isolated skyrmion 
falls off with distance only as r^ 1 . We argue in ap- 
pendix |^ that long range ferromagnetic order is indeed 
preserved, and that the internal degrees of freedom of 
dilute skyrmions and antiskyrmions develop the correla- 
tions necessary to ensure this. 

In contrast, at strong disorder (A ^> J), the charge 
density provides almost perfect local screening of the dis- 
ordered potential, so that 

Sp(r) ~ -V(t)/U (A/J»l). (11) 

Corrections to perfect Thomas-Fermi screening arise at 
length scales larger than Ln, where exchange becomes 
important. The effect of exchange is to force screen- 
ing charges to be quantised, since an unquantised charge 
costs divergent exchange energy in the thermodynamic 
limit. We can summarise the effect of exchange by divid- 
ing the system into regions of area L^, finding for each 
such area the integral 

Q=~ I ^ (12) 
Jl 2 h Uq 

and adjusting the total screening charge within every re- 
gion to the integer value closest to Q. We argue that 
these integers are predominantly zero in a ferromagnetic 
phase, and predominantly non-zero in a phase without 
ferromagnetic order. To see this, consider a well-ordered 
ferromagnetic phase, in which S(r) has small spatial vari- 
ations around a global direction of magnetisation. In this 
case, the net topological charge in any region has magni- 
tude much less than one. Conversely, in a phase without 
such order, unit topological charge will typically accumu- 
late over a region of linear size given by the ferromagnetic 
correlation length. 
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This picture leads us to identify the phase boundary of 
the ferromagnet as the point at which (Q 2 ) 1 / 2 ~ 1. To 
estimate (Q 2 ), note from Eq. jl^) that each correlation 
area of size A 2 contributes to Q a charge of magnitude 
X 2 A/Uo and random sign. Over the Hartree area, there 
are L 2 H /\ 2 such contributions, so that 

(Q 2 ) 1/2 ~ (X 2 A/Uo)(L H /X) = AL H A/{/ . (13) 

The phase boundary is therefore given by 

A c ~ U /(XL H ) = J(L H /A). (14) 

For A > A c , the ground state is strongly disordered and 
has no ferromagnetic order. Within our classical descrip- 
tion, the spins are frozen at zero temperature. We can 
therefore identify this phase as a quantum Hall spin glass. 

B. Away from integer filling 

Consider next the ground state away from v = 1. We 
examine behaviour at fixed A < A c , MS (X function of 
the average charge density (Sp) oc v — 1. For (dp) = 
and A < A c , the system has a net magnetisation, but 
the introduction of charge, in the form of skyrmions (or 
anti-skyrmions) disrupts ferromagnetic spin alignment. 
The size R of an isolated skyrmion in a clean system 
is divergent if Hartree repulsion is not balanced by Zee- 
man energy. However, the presence of an impurity poten- 
tial establishes an optimal size, because it is energetically 
favourable for a skyrmion to locate its charge distribution 
in randomly occurring potential wells. The typical depth 
of such a well of radius R 3> A, averaged over its area, is 
— A(A/i2). (The case A > R is treated in Ref.[lO|) The 
Hartree energy is ~ Uq/R 2 . The total energy is hence 
minimised if 

R~U /(AA) = L H (A C /A). (15) 

Note that this value of R exceeds Lh for A < A c . Be- 
cause of this, we expect that exchange energy dominates 
over the potential energy and that the skyrmion will not 
be strongly distorted by the random potential. In other 
words, screening as discussed in the previous subsection 
does not alter the present argument. 

We expect that ferromagnetic order will persist with 
increasing charge density until such skyrmions overlap. 
The phase boundary hence lies at 

(^) C = L H 2 (A/A C ) 2 . (16) 

In summary, we have used simple arguments to ob- 
tain the phase boundaries between the ferromagnetic and 
the spin glass as a function of disorder and charge den- 
sity. The results, Eqs. ( |l4] ) and (|l6|), are summarised in 
the schematic phase diagram shown in Fig.|l|. We next 
present results from Monte Carlo simulations which sup- 
port these predictions. 
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FIG. 1: Schematic diagram for Lh/A S> 1. 

The discussion we have presented is for a model 
with a short-range Hartree interaction and for the limit 
is < A < La. The central consequence of using a 
Coulomb form, U(r) = e 2 /47reor, in place of a short- 
range Hartree interaction is to change the length scale 
derived by comparing Hartree with exchange energies, 
from L H = (U /J) 1/2 to L H = e 2 /47re J = 16(2tt) 1 / 2 / s . 
Since in this case Lh is not parametrically larger than I b , 
scaling arguments of the type presented in this section 
are not justified. In this case, in place of three distinct 
regimes of disorder strength (A « J, J « A < A c and 
A c <C A) we expect only two, with A c ~ J. 

IV. MONTE CARLO SIMULATIONS 

In order to test the theoretical results derived in the 
previous section, we use Monte Carlo methods, in com- 
bination with simulated annealing, to study a classical 
Heisenberg model with quenched disorder. In this sec- 
tion, we outline our methodology. Similar techniques 
have been applied recently to systems without quenched 
disorder in Ref.[l6|. We present our results for the phys- 
ical response functions of the model in sections [v] and 

v| 

We treat a lattice version of Eq. (||) with NxN classical 
spins Si of unit length on a square lattice, taking peri- 
odic boundary conditions. Spins have nearest-neighbour 
ferromagnetic interactions of strength J. In addition, the 
electrical charge qk on each plaquette k is calculated from 
the area on the spin sphere covered by the spins on the 
corners of that plaquette. This charge density has a local 
repulsive Hartree interaction of strength Uq, and is also 
subject to a uniformly distributed random background 
potential e/- S [—A, A]. Since we choose this potential 
independently for each plaquette, the correlation length 
A is set by the lattice spacing. 

The Hamiltonian of the lattice model is 

(ij) k 

+7(Qo-$>fc) 2 , (17) 
k 
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where we have introduced a Lagrange multiplier 7 to bias 
the system towards a predefined number Q of charge 
quanta. 

To obtain the ground state spin configuration for a 
given disorder realisation, we start from a random ini- 
tial state and anneal using Monte Carlo dynamics and 
the Metropolis algorithm. After some experimentation, 
we found the following three-stage protocol to be effec- 
tive. In the first stage, the temperature is reduced lin- 
early in time from a high temperature To (several times 
J) to Tq/10, using 3.5 x 10 5 Monte Carlo steps per spin 
(MCS). In the second stage, the temperature is reduced 
from To/10 to 0, using 5 x 10 5 MCS. For both these 
stages, the attempted spin update is an isotropically dis- 
tributed reorientation. In the third stage, the system is 
quenched for a further 5 x 10 4 MCS, using as the at- 
tempted spin update only small angle reorientations in 
order to improve the acceptance rate. 

We have checked that this algorithm finds the ground 
state reliably for a weakly disordered system with overall 
charge neutrality, by doing repeated runs for a given re- 
alisation of background potential and using local charge 
and energy densities to identify states that differ only 
by a global spin rotation. In strongly disordered sys- 
tems and those with non-zero average charge density, re- 
peated applications of the algorithm do not reproduce 
the same state to high precision. Instead, a number of 
low-lying states are obtained, having a small spread (less 
than 10%) in their energies and other observables such as 
their magnetisation. We perform at least five indepen- 
dent simulations for every disorder realisation and pick 
from the states obtained the one with the lowest energy. 
For a chosen disorder strength A, we generate three or 
more realisations of the disorder potential and average 
our results over these realisations. We find that proper- 
ties such as the ground state energy, the magnetisation 
and various response functions show only small fluctua- 
tions between different disorder realisations for systems 
of size 32 x 32 spins or larger. Simulations of very much 
bigger systems are ruled out by constraints on computing 
time, and we choose N = 40 as the system size for most 
of our simulations. 

The Lagrange multiplier 7 is chosen sufficiently large 
that a majority of the simulation runs yield a ground 
state with the desired total charge Qq. We do not detect 
any dependence of our results on the particular value of 
7- 



A. Magnetisation 

We calculate the site-averaged magnetisation 

M=\{S(r))\=N- 2 \J2Si\. (18) 

i 

Results are shown for a system of 40x40 spins with 
Uo = 8nJ in Fig.||(a). Spins are fully aligned for 
A/4ttJ < 1, because our bounded disorder distribution 
then lies entirely below the threshold, Eq. Increas- 
ing A beyond 47r J, we obtain a partially polarised ferro- 
magnet. In fact, the magnetisation remains close to its 
saturated value until A/47rJ ~ 1.8, when it starts to drop 
appreciably. Increasing the disorder strength further, we 
reach a regime of strong disorder at A/47rJ ~ 3, in which 
there is only a small magnetisation, of the magnitude ex- 
pected for a spin-glass state of a finite-size system. 

These results are consistent with the existence of a 
phase transition from the collinear, ordered phase at 
weak A to a spin-glass phase at strong disorder. From 
the magnetisation curve we pick A c /4irJ ~ 2.5 as our 
estimate of the critical disorder strength. 
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FIG. 2: Ground state properties for a system with Uo = 8irJ 
as a function of disorder strength: (a) magnetisation M as 
a fraction of the saturated moment, (b) susceptibility \, (c) 
spin stiffness p a in units of J. 



V. MAGNETIC RESPONSE 

Having outlined our simulation techniques, we next 
present our results. We start with the magnetic prop- 
erties of the system. We show that they indicate a tran- 
sition from the ferromagnet to the spin glass with increas- 
ing disorder strength, as expected from our discussions 
in section III. 



B. Susceptibility 

We also calculate the uniform susceptibility \ from the 
response of the ground state to a weak Zeeman field. Our 
procedure is as follows. First we obtain a ground state 
without a Zeeman field, using the protocol described in 
Sec. IV. Then we apply a weak Zeeman field h in the 



direction of the residual magnetisation, by adding to the 
Hamiltonian the perturbation SH = —h ■ M. We find 
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the ground state in the presence of this perturbation by 
repeating the third stage of the quenching protocol. Tak- 
ing care to check that our measurements remain within 
the linear response regime, we extract the susceptibility 
from the change in the magnetisation 5M , using 



X = \SM\/\h\. 



(19) 



Results for Uq — 8ttJ are shown in Fig. |2j(b) . The suscep- 
tibility as a function of disorder strength has a large peak 
at A/AttJ = 2.5. We interpret this as a second indication 
of a phase transition from the ferromagnetic phase to a 
disordered phase. 



C. Spin Stiffness 

We now turn to the spin stiffness p s which measures 
the rigidity of the spin configuration. It is obtained 
by calculating the energy cost of small amplitude, long- 
wavelength spin twists in the ground state. More specif- 
ically, let the columns of sites in the N x N lattice 
be labelled by integers 1,2, ...N. Then, starting with 
a ground state obtained as in Sec. IV, we construct a 
twisted state by rotating all spins on column N/2 of the 
system through a small angle 8 about an axis e. Using 
the third stage of the quenching protocol, we then relax 
all spins in this twisted state except those on columns 
1 and N/2, which are held fixed. From the difference in 
energy AE between the initial and final states, we obtain 
the spin stiffness for rotations about the axis e, using 



p s = AE/26 2 



(20) 



Repeating this for different axes of rotation, we calcu- 
late the full 3x3 symmetric tensor for the spin stiffness, 
p"^ '. As expected, in the ferromagnetically ordered phase 
one of the principal axes of this tensor lies to a good ap- 
proximation along the magnetisation direction, and it is 
convenient in these calculations to choose rotation axes 
e in directions parallel and perpendicular to M. 

The eigenvalues of the spin stiffness tensor are shown 
as a function of disorder strength in Fig.|^(c). For the 
fully-polarised ferromagnet, rotations about the magneti- 
sation direction do not alter the spin configuration, and 
so one eigenvalue is zero while the other two are degener- 
ate, taking the value p s = J. For the partially-polarised 
ferromagnet, all three eigenvalues are non-zero, with two 
remaining degenerate. The spin-glass state, however, 
has no special spin direction and this magnetic isotropy 
means that that all three eigenvalues are approximately 
degenerate. The stiffness is reduced in value from J in 
the fully polarised ferromagnet to approximately J/2 in 
the spin glass. (A variational estimate for the stiffness in 
the spin-glass phase is 2 J/ 3. See appendix |a|) The mag- 
netically isotropic phase is observed for A /Air J > 2.5, 
yielding the same estimate of A c as our magnetisation 
and susceptibility data. 



D. Spin Correlation Length 

The behaviour of spin correlations provides a further 
way of characterising ground states. In particular, we 
consider the correlation function 



C(r) 



(r) 



(21) 



where the sum runs over all N r spin pairs of separation 
r. We extract the spin correlation length £ by fitting to 
the form 



C'(r) = M 2 + (1 - M 2 ) exp(-r/£) 



(22) 



The behaviour of the correlation length as a function 
of the disorder strength in the spin-glass phase is shown 
in Fig.|3[ From the Harris criterionEZI, we expect it to 
diverge as £ ~ (A — A c )~ v with v > 2/d = 1 as the ferro- 
magnetic phase boundary is approached. Our results are 
consistent with a divergence at A/AnJ = 2.5, although 
it appears that our they are affected by finite-size effects 
for £ for A/4irJ < 2.7. Perhaps because of these finite- 
size effects, this fit gives a low value for the exponent: 
v = 0.7. Attempts at a similar fit in the ferromagnetic 
phase are unsuccessful, and indeed the arguments of Ap- 
pendix |c] suggest that correlations in this case may decay 
with a power law. 
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FIG. 3: Correlation length (in units of the lattice constant) 
as a function of the disorder strength. Open symbols: our 
calculation, solid line: fit to a power law. 



E. Phase Diagram 

Our Monte Carlo results outlined above allow us to dis- 
tinguish between a ferromagnetic phase and a quantum 
Hall spin glass. For a system with Uq = 8irJ, we conclude 
that a phase transition from a collinear ferromagnet to a 
spin-disordered phase occurs at A/AttJ — 2.5. Repeat- 
ing these calculations for different disorder strength and 
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charge density, we can map out a phase diagram for the 
ground state of the system. Results for U — Air J are 
shown in Fig. El 




A / 471 J 



constant e(q) or by the compressibility «((/), related via 

K(q) = q 2 e(q)e /e 2 , (23) 

where e is the electron charge and eo is the permittiv- 
ity of the medium. More precisely, we apply a periodic 
modulation, V(r) — > V(r) + Vx cos(q ■ r), to the potential 
in Eq. (Q). As a result, the ground-state electron den- 
sity changes according to 5p(v) — ► Sp(r) + 6pi(r). Since 
the system is disordered, Spi(r) contains many Fourier 
components, but after averaging the linear response is 

(5pi(r))=-Ki/t(g)cos(qT), (24) 

which constitutes our definition of n(q). 

Our numerical results for n(q) are displayed in Fig.0, 
where we study systems deep in the spin-glass phase, 
setting A/ An J = 3 and (8p(r)) = 0. We compare be- 
haviours at Uo/AttJ = 1 and at Uo/AirJ — 2, in each case 
combining data from lattices of size 40 2 and 56 2 in order 
to maximise wavevector resolution. 



FIG. 4: The phase diagram for Uo = An J (Lh/A = VAn), as 
a function of disorder strength and charge density. QHFM: 
quantum Hall ferromagnet, QHSG: quantum Hall spin glass. 
Inset: Results from scaling arguments for Lh/A 2> 1. 

Comparing this with the phase diagram (Fig. |l|) pre- 
dicted from our heuristic arguments, we see that the two 
are qualitatively very similar, even though the simula- 
tions are carried out for (L#/A) = ^/An while our scal- 
ing picture applies in the limit (L#/A) ^> 1. An idea 
of the dependence on Lh/X is given by contrasting re- 
sults at the two values of Uo studied. The critical disor- 
der A c at which the neutral system loses ferromagnetic 
order is A c /AnJ ~ 2.5 for Uo = 8irJ, and reduces to 
A c /AirJ ~ 2.2 for Uo = AirJ, in qualitative a gree ment 
with the scaling behaviour we expect from Eq. (|l4J). 

It is interesting to note that there is a density range 
over which disorder may stabilise the ferromagnet, by 
limiting the size of the nucleated spin textures. This 
range is however very narrow: \(Sp) \ < ICR 2 in units of 
charge per plaquette. 

We mention in passing that we have not searched ex- 
tensively for a skyrmion crystal, expected at a finite 
charge density in the weak disorder limit but presumably 
unstable for A ^ 0. 



VI. DIELECTRIC RESPONSE 

We have so far discussed the phase diagram of the sys- 
tem in terms of its magnetic correlations and response. 
We next study its charge response. 

We examine first the dielectric response of the par- 
tially polarised ferromagnet and spin glass, characterised 
at zero frequency by the wavevector-dependent dielectric 
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FIG. 5: Compressibility n(q) as a function of wavevector q 
for systems with J = 1 and Uo ~ 4n and 8-7T. (A/AttJ — 3, 
(dp) = 0.) 

We find that, at small q, n(q) is independent of Uo 
and quadratic in q. At large q, it is independent of q, 
varying roughly as U^ 1 . We can understand these results 
for n(q) using the approach we employed to discuss the 
phase diagram. The Hartree length in again plays an 
important role. 

For q 3> L^ 1 , exchange may be neglected and we see 
from Eq. |J|) that 

K (q) ~ U -' (qL H » 1) , (25) 

in agreement with the q-independent value obtained for 
n(q) at large wavevectors from our simulations, and with 
the Uq dependence of these values. 

Exchange becomes important at longer wavelengths 
q ijj 1 . To estimate (6pi(r)) for small q, we suppose 
that it arises primarily from spin rotations which have 
amplitude 69 and wavevector ~ q. To be specific, con- 
sider for the periodic perturbation V\ cos(qx) a region of 
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FIG. 6: Mapping a line of spins on the lattice to a trajectory 
in spin space. 

size 2L x L with L = Tt/q. The perturbation causes a 
net movement of charge SQi in the x direction, from one 
side of this region to the other. We choose to focus on 
those spins lying on the line of length L which is paral- 
lel to the y axis and which divides the region into two 
equal halves of size L x L. (See Fig. §) We suppose, in 
a variational spirit, that the charge movement induced 
by the perturbation involves simply a rigid rotation of 
these spins, through an angle 69. With this assumption, 
let us estimate the charge transfer generated by such a 
rotation, and the exchange energy that it costs. To do 
so, it is useful to regard the spin configuration as a map 
of the line we have defined onto a trajectory on the spin 
sphere. For a disordered spin configuration, this trajec- 
tory resembles a random walk. The end points of the 
walk are randomly placed on the spin sphere if L is large 
compared to the spin correlation length £, and the dis- 
tance D between them on the surface of the spin sphere is 
therefore typically 0(1). Under a rigid rotation of spins 
on the line, the trajectory is displaced rigidly around the 
spin sphere. The charge transferred across the real-space 
line is proportional to the area swept out on the spin 
sphere by the trajectory during this displacement. We 
therefore arrive at the estimate 

6Qi~D59~ 59. (26) 

Choosing the phase and axis of rotation appropriately 
in each such region, the change in potential energy den- 
sity arising from the rotation is — ViSQi/L 2 = —V\q 2 59 
while the associated change in exchange energy density 
is Jq 2 59 2 . Choosing 59 to minimise the total energy, we 
find |<5pi(r)| ~ V\q 2 / J and hence 

K (q)~q 2 /J = U 1 (qL }i ) 2 [qL H « 1) (27) 

This conclusion is again consistent with our numerical 
results. 

Alternatively, we can arrive at this form for the com- 
pressibility from scaling considerations. In general, we 
may expect that the compressibility is described by the 
scaling form 

K{q) = U^f{qL H ) , (28) 

where f(x) approaches a constant at large x. Our central 
hypothesis is that exchange dominates at small wavevec- 
tors, and so n(q) should be independent of Uq as q — > 0. 



This implies that f(x) ~ x 2 for x <§; 1, leading to the 
form given in Eq. (E7|) for the compressibility at small q. 

To summarise, the system has a metallic response to a 
perturbing potential at large wavevectors, with n(q) inde- 
pendent of q. However, at small wavevectors, it behaves 
like an insulator, with e(q) independent of q. 

We note that, for an infinite system at exactly q = 
and zero frequency, /c(0) should be proportional to the 
thermodynamic density of states. So, except in the fully- 
polarised ferromagnet, we expect that k(0) remains fi- 
nite as we take the temperature to zero. To reconcile 
this general expectation with our results, one must re- 
member that for finite-size systems considered here at 
zero temperature, the discreteness of charge implies that 
almost all disorder realisations have no response to an in- 
finitesimal change in the chemical potential, while there 
is a divergent response from those realisations for which 
ground states from two different charge sectors are de- 
generate. 



VII. OPTICAL CONDUCTIVITY 

Finally, we consider the optical conductivity a(u) at 
frequency u). Within our treatment, spin waves are the 
only excitations that contribute. Since they arc not 
topological excitations and so do not carry net electri- 
cal charge, the dissipative conductivity vanishes in the 
low- frequency limit. Spin waves do, however, give rise to 
local charge fluctuations and a fluctuating electric dipole 
moment which couples to an oscillating external electric 
field, generating dissipation at finite frequency. 

For.— a fully-polarised quantum Hall ferromagnet, 
Greenlia has shown that the spin-wave contribution to 
optical conductivity is very small. In the non-collinear 
quantum Hall spin glass, the contribution may be larger 
due to the presence of a finite charge density in the disor- 
dered ground state. Furthermore, the low-energy dynam- 
ics of the collinear quantum Hall ferromagnet is qualita- 
tively different from that of a quantum Hall spin glass, 
since while spin waves in a collinear background have 
a quadratic energy dispersion, those-jn a non-collincar 
background have a linear dispersionlia, uj = cq at small 
q, with velocity 

c=( Ps /x) 1/2 - (29) 

Of the three polarisation modes, one is expected to re- 
main gapless even in the presence of a Zeeman coupling. 

We now calculate the spin- wave contribution to the op- 
tical conductivity for a ground state with non-collincar 
spins, and then estimate its magnitude in realistic sys- 
tems. 

The rotation of spins from their ground-state orienta- 
tion in the presence of a spin wave may be parameterised 
by a vector p(r), with at first order 

S^S-p{r)xS. (30) 
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Using Eq. (||) , we show in appendix [A| that this induces 
a change in the electron density 

SpW^Sptf-^dip-djS. (31) 

This linear dependence on Vp should be contrasted with 
spin waves in a collinear ferromagnet, for which the den- 
sity change is 0(|Vp| 2 ). 

To estimate the contribution to the conductivity, we 
note that a spin wave will couple to an external electric 
field through the net dipole moment P that it induces. 
The coupling appears in the Hamiltonian in the form 

SH = -P x £ x (32) 

for an electric field £ x in the x-direction. From Fermi's 
golden rule, the power H(ui) absorbed from an oscillating 
electric field of frequency to is 

= -^£ 2 \(f\P x \0)\ 2 hu;g(cj) (33) 
n 

where |0) is the ground state and |/) is a state with a 
single spin wave excited, of excitation energy %u>, and 
g(u>) is the spin wave density of states in frequency. From 
Eq. (|29|), we have 

9(^) = ^L x L y , (34) 

including three polarisations for a system with linear di- 
mensions L x and L y . 

We estimate for the matrix element appearing in 
Eq. (|3^) by considering the dipole moment induced by 
a spin wave. We start our discussion using the form 
taken by a spin wave mode in the absence of ground- 
state disorder, p(r) = ap cos(q • r), where p is a unit 
vector defining the axis of spin rotations and a specifies 
the amplitude. The charge density induced by an excita- 
tion of this type in a disordered ground state with spin 
correlation length £ has, from Eq. (|3l|), magnitude aq/£- 
We expect this to fluctuate with a random sign over the 
length scale £. The spin wave therefore induces electric 
dipoles of magnitude e£,(q£a) in each correlation area £ 2 . 
Averaging over aijX L y system, we find 

<P 2 > ~ (e£V) 2 ^ ~ eV (jff^Ly. (35) 

where we have substituted q = uj/c. 

It remains to determine the amplitude a for a single 
quantum excitation. We show in appendix [b] that 

(a 2 ) = h/ X LuL x L y . (36) 

Combining factors and dropping numerical coefficients, 
the absorbed power is 

nH-^^. (37) 

Ps 



This is an ohmic response, H/L x L y — er£ 2 /2, with con- 
ductivity 

. . e 2 / uj£ \ 2 fioj ,„„ N 

w~t{t) Ts (38) 

Let us estimate the magnitude of this spin wave con- 
ductivity. In the spin glass phase, we use our numerical 
results to estimate p s ~ J/2. Also, the numerical results 
for the lattice shows that xiattice of the order of unity in 
units of 1%/J- The continuum magnetisation is related 
to the lattice spins by: fa <-> Sfi/l 2 B . This means that 
the continuum susceptibility is x — h 2 / Jl%- Combining 
these factors, we have 

Taking £ = VMb and J = 4-8 kelvin (correcting 
our earlier valudj), we find for a frequency of 1 GHz, 
<j(gj) w p,0 _3 -10 -4 ) e 2 /h. Unfortunately, variable-range 
hoppingE] seems likely to mask this contribution to a(ui). 

VIII. SUMMARY AND DISCUSSION 

We have investigated the competition between ex- 
change interactions and disorder in quantum Hall fer- 
romagnets at or near integer filling and at zero tempera- 
ture. Our approach is tailored to the limit of a smoothly- 
varying impurity potential. 

We find that the ferromagnetic state is destroyed by 
strong disorder through the creation of skyrmion/anti- 
skyrmion pairs, or by a finite density of either skyrmions 
or anti-skyrmions at filling factors sufficiently far from 
v = 1. This behaviour, anticipated from simple scaling 
arguments, is confirmed in numerical studies of ground 
state spin configurations, obtained from slow Monte 
Carlo cooling of an initially random high temperature 
phase. The disordered phase is identified as a quantum 
Hall spin glass by the absence of long-range order, the 
presence of non-vanishing local magnetic moments and a 
finite spin stiffness. 

The quantum Hall spin glass has a zero-frequency di- 
electric response which interpolates between that of an 
insulator at small wavevectors, and that of a metal at 
large wavevectors. It supports gapless spin wave modes 
that couple to electric fields through a finite dipole mo- 
ment and contribute to the optical conductivity of the 
system. 

Possible experimental signatures of the phenomena 
we have discussed follow both from the behaviour of 
the magnetisation and from the nature of excitations. 
Measurements of fthfi-Knight shift in nuclear magnetic 
resonance (NMR)EifEj provide information on the dis- 
tribution, sampled in space, of the spin polarisation 
component parallel to the applied magnetic .fifild, while 
polarisation-resolved absorption spectroscopyBcJ can be 
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used to determine the average spin polarisation. Impu- 
rity effects are most characteristic at v = 1, where they 
result in a reduced spin polarisation, reaching zero in the 
spin glass at zero Zeeman coupling, and a finite width 
to the polarisation distribution. We note that spin po- 
larisation which remains unsaturated even at v = 1 is 
found in absorption spectroscopyj23'E-J and that a broad 
Knight-shift distribution is measured in low-temperature 
NMR.)E3 In addition, the existence of excitations at en- 
ergies lower than the Zeeman gap is characteristic of the 
partially polarised ferromagnet, and has been advanced 
as an explanation for a coupling observed in NMR exper- 
iments between padio-frequency magnetic fields and the 
electron system.ES On the other hand, at filling factors 
away from v = 1, neither local probes nor the nature of 
excitations distinguish sharply between a skyrme crystal 
in a clean system and the partially polarised ferromagnet 
or spin glass induced by disorder. 

We are grateful for discussions with N. R. Cooper and 
S. L. Sondhi. The work was supported in part by the EP- 
SRC under Grant GR/J78327 (JTC), and by the Royal 
Society (DKKL). 



APPENDIX A: SPIN WAVES AND CHARGE 
FLUCTUATIONS 

Consider a small rotation R(p) of the spins around p, 
so that spin directions transform according to S — ► RS, 
with 



R 



and 



[P] a6 



e abc p c 



(Al) 



To second order, SS ~ -p x S + p x (p x S)/2. The 
charge density deviation from v = 1 is given by Spo = 
e-ijS ■ (diS x djS)/8ir. In the presence of an additional 
spin rotation, it becomes 



so that 



S Pl = ^RS- 



d l {RS) x dj(RS) 



(A2) 



S/n = -ji-dip- d,S + ^dip- 3, {p xS) + 0( P 3 ) , (A3) 

47T 87T 

obtained using (R _1 <3jR)5* ~ (d t P - [P,diP]/2)S = Sx 
dip + (dift x p) x 5/2. We can drop higher-order terms if 
|p| < 1 and |Vp| < |V5|, so that Spi/8p < 1. 

To first order in p, the component of p parallel to S 
does not affect the charge density. At this order, the 
charge fluctuation Spi can be written as 

8 Pl = -e ij d i p-d j §. (A4) 

From continuity, Spi + diJi = 0, the current density is 

Ji^tijP-djS + divergence-free part (A5) 

Note that we have only identified the transport current, 
and not any circulating current in the bulk. 



Let us now consider the energy cost of spin rotations 
for the Hamiltonian in Eq. (Q) . The change in the ex- 
change energy density Hj is 



SHi 



J 2 



Jd l p-{d l Sx S)+8H ( f ] , 

|Vp1 2 - 2(Vp • S) 2 + {dip- S)8 % {p- S) 



-(p-SXdtf-diS) 



(A6) 



The change in potential energy density is 



SH p = V n (r)S Pl (r) 

where 

V h (t)=V{t) 



^i(r)[/(r-r')^i(r')d 2 r' 
(A7) 



ft+ U{r-r')5p (r')d 2 r' (A8) 



is the local Hartree potential and p is the chemical po- 
tential. 

An equation satisfied by the spin configuration in the 
ground state is obtained from requiring that the total 
energy is unaffected by the rotation p to first order. This 
means that we have to balance the first-order terms in 
the expressions for SHj and SH p , giving 



0, 



Jd t S x S -cijVni^djS 



= 



(A9) 



for the spins in the ground state. This non-linear equa- 
tion demonstrates the competition between ferromag- 
netic exchange and Thomas-Fermi screening in a disor- 
dered QHFM. 

(2) 

From the second-order contributions SHj in the ex- 
change energy SHj , we can make a variational estimate of 
the spin stiffness in the spin-disordered phase. Suppose 
that p and S are uncorrelated. Then the second-order 
terms average to 

1 



((Vjf • sy) = {{d^ ■ s)di& ■ s)) = -ivpi 

((p-S)(dip-diS)) = 0. 



(A10) 



We can define a disorder-averaged spin stiffness from 

12 , 



(SH ( j 2) ) =p s |Vp1 2 /2 so that 



2 J/3. 



(All) 



APPENDIX B: LONG WAVELENGTH SPIN 
WAVES 

In this apoendix, we review the results of Halperin 
and Saslpailj for hydrodynamic excitations, and of 
GinzburgEEl for elementary excitations, and adapt these 
for our purposes. 

Both theories deal with long- wavelength, low-energy 
excitations of a disordered spin system. The microscopic 
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magnetisation density and energy density can be coarse- 
grained over areas A to mean values; in the case of the 
magnetisation one has 



1 



ieA 



(Bl) 



Within a hydrodynamic theory, the low-energy dynam- 
ics of the system is determined entirely these coarse- 
grained quantities, which obey conservation laws and are 
assumed to fluctuate slowly on the timescale set by local 
relaxation rates. 

Hydrodynamic spin fluctuations have equations of mo- 
tion which may be written in terms of the rotation angle 
p, as used in in Eq. (|30|), and the local magnetisation, 
with the form 



dp a _! 
dm a 2 



(B2) 
(B3) 



where p s is the spin stiffness and y is the uniform mag- 
netic susceptibility. Using the same variables, the free 
energy of the system is 

AF(rn,p) = ipJd 2 r(±- m i+p s \V Pa \^ 

= / d 2 r( X pl+Ps\Vp a \ 2 ) . (B4) 

a=l J 

A long-wavelength treatment of elementary excitations 
leads to the same equations of motion and an equivalent 
expression for the energy of the system, in which p s and y 
are ground-state quantities. We have argued in section |\] 
that both are finite at zero temperature in the spin-glass 
phase. These equations are reminiscent of the dynamics 
of the spin waves in a Heisenberg antiferromagnet. They 
lead to spin waves with the dispersion relation 



to 



cq; 



(Ps/X) 



1/2 



(B5) 



Considering a single spin-wave mode, with 

p = ap cos(q ■ r - cut) (|p | = l), (B6) 

we obtain the amplitude a for one quantum by setting 
the energy to huj, arriving at Eq. (pq). 



APPENDIX C: DILUTE SKYRMIONS 

In this appendix, we study in more detail the interac- 
tions between a dilute set of skyrmions at weak disorder. 
We show how the relative orientations of internal degrees 
of freedom of the skyrmions are determined by the ener- 
getics of the system. 

For a clean system without Hartree interactions, the 
Hamiltonian reduces to the 0(3) non- linear sigma model. 



The skyrmions are of the Belavin-Polykov type. They do 
not interact. Hartree interactions alone will lead to a di- 
vergent skyrmion size in the absence of a Zeeman field. 
This is prevented in the presence of a disordered poten- 
tial, as pointed out in section HI. However, the pres- 
ence of an inhomogeneous potential also means that the 
skyrmions interact. We discuss here whether a physical 
picture of isolated skyrmions at weak disorder is justified. 
In particular, since the spin deviation due to a skyrmion 
falls off as 1/r with distance r, one should ask whether 
the ferromagnetic polarisation is strongly affected by a 
collection of dilute skyrmions. 

For the purposes of this appendix, we choose a sim- 
ple disorder distribution in which pinned skyrmions and 
anti-skyrmions are nucleated by isolated wells and barri- 
ers of circular shape, positioned randomly on the plane. 
The depth (or height) of these potentials exceeds the 
threshold, Eq. (^), so that they nucleate skyrmions (or 
anti-skyrmions). This model disorder distribution has 
the advantage that we can use the Belavin-Polyakov so- 
lutions of the non-linear sigma model as a starting point 
for our analysis. 

Consider first the case in which the chemical poten- 
tial is sufficiently high that there are no anti-skyrmions. 
A set of Belavin-Polyakov skyrmions is described by an 
analytic function w(z) of the position z = x + iy in the 
plane 



w( Z )=n- 



(Cl) 



where n is the number of skyrmions. The function w(z) 
parametrises the spin configuration via 

S x +iS y = 2w/(l + H 2 ) 

s z = (i-M 2 )/(i + M 2 ) (C2) 

with the boundary condition S x — 1 at infinity. For a 
single skyrmion (n — 1), the density profile is 



5p n =i(z) 



\d z w\ 



d?/4rr 



tt(1 + M 2 ) 2 (\ z - Zi \2 + dl/4y 



(C3) 



where Zi — (dj + 6j)/2 can be interpreted as the position 
of the skyrmion and di = |a$ — bi\ as its diameter. There 
remains an internal phase degree of freedom (f>i, defined 
from (a, — hi) = diexp{i<pi). In multiple skyrmions con- 
figurations the relative phases are important. 

In a dilute glass, skyrmions are typically separated by 
distances large compared to their size, and the position 
and diameter of each skyrmion is determined separately 
by the balance of Hartree and potential energy in the 
vicinity of the potential well on which it is centered. In 
the region far from any well, the Belavin-Polyakov form 
is a good description for the spin configuration, because 
it minimises exchange energy. The parameters Zi and di 
are fixed by the spin configuration within the ith well, 
while the values of 6; remain to be determined. The 
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density profile in this region is 
Sp'(z) = 



7T(1 + M 2 ) 2 



E 



(z - a,i)(z - bi) 



(C4) 



where the prime indicates the restricted area in which 
this form applies. The interference terms in the squared 
sum make a contribution to the total charge in the re- 
gion far from the wells which amounts to a pairwise cou- 
pling between the internal degrees of freedom of different 
skyrmions. The integrated density is 



I = 



7T (1 + H 2 ) 2 



(z* -a*)(z* -b*)(z-a 3 ){z-b 3 ) 



(C5) 



In the dilute limit, we can take a; 



Zi. We can also 



approximate \wy 
I ~ Re 



1 away from the skyrmions. Then 



1 d 2 



z \ - 



4tt ^ (z* 



1 x didj cos(4>i — 4>j) 
2^ \zi- Zj \* 

1>J 1 J 



(C6) 



This contribution to the charge in the region between 
potential wells adds to the energy of the system, since 
the electrostatic potential here does not exceed thresh- 
old [Eq. (g)]. It represents charge density that is removed 
from the cores of skyrmions, as a result of overlap with 
the tails of other distant skyrmions. We see that in order 
to minimise the energy we should choose the skyrmion 
phases <pi so as to minimise /. This is done by arranging 
4>i so that the variable die 1 ®* sums locally to small val- 
ues. With such correlations, we expect that long-range 
ferromagnetic order survives in the presence of dilute 
skyrmions. 

A similar treatment of a mixed system containing both 
skyrmions and antiskyrmions leads to interactions of 
the same kind amongst the skyrmions, and separately 
amongst the antiskyrmions, but without coupling be- 
tween the two species at leading order in inverse density. 
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